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Estimated  UV  clutter  levels  at  10-100  meter  sensor  pixel  resolution 
extrapolated  from  recent  Polar  Bear  measurements 

M.  Wohlers,  R.  Huguenin,  and  M.  Weinberg 

Aerodyne  Research,  Inc. 

Billerica,  MA  01821 

and 

R.  Huffman,  R.  Eastes,  and  F.  Delgreco 
Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731 


ABSTRACT 

There  is  continuing  need  for  information  about  the  earth  background  clutter  at 
ultraviolet  wavelengths.  This  paper  describes  the  methodology  and  the  results 
obtained  at  1304  A  wavelength  from  an  analysis  of  the  AFGL  Polar  Bear  experiment. 
The  basic  measurement  equipment  provided  data  of  a  spatial  resolution  of  20  km  over 
a  large  portion  of  the  earth.  The  instrumentation  also  provided  sampled  outputs  as 
the  footprint  scanned  along  the  measurement  track.  The  combination  of  the  fine 
scanning  and  large  area  coverage  provided  opportunity  for  a  spatial  power  spectral 
analysis  that  in  turn  provided  a  means  for  extrapolation  to  finer-spatial  scale;  a 
companion  paper  discusses  the  physical  basis  for  this  extrapolation. 
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There  is  continuing  need  for  information  about  the  earth  background  at 
ultraviolet  wavelengths.  There  is  currently  little  or  no  experimental  data  from 
which  we  can  infer  the  spatial  variability  of  ultraviolet  earth  backgrounds  at 
spatial  scales  below  1  km.  While  it  is  evident  that  new  flight  experiments  will  be 
required  to  obtain  this  needed  background  data,  design  of  the  new  experiments  will 
require  some  knowledge  of  the  background  radiance  properties  at  the  scales  that 
will  be  measured.  Clutter  irradiances  at  the  desired  scales  are  needed  to  insure 
that  the  sensor  design  includes  an  adequate  optical  system,  detection  response,  and 
measurement  scheme. 


Recognizing  this  apparent  dilemma,  we  examined  the  feasibility  of  extracting 
finer  scale  spatial  clutter  characteristics  from  existing  UV  backgrounds  data.  UV 
backgrounds  measurement  programs  were  examined,  and  a  significant  opportunity  to 
extract  the  desired  information  from  Polar  BEAR/AIRS  imager  data  was  discovered. 
Because  of  the  unique  measurement  sequence  employed  in  the  experiment,  the 
potential  existed  for  extracting  clutter  information  at  spatial  scales  as  fine  as 
45  meters.  The  significance  of  this  possibility  is  increased  by  the  fact  that  the 
Polar  Bear/AIRS  Imager  is  currently  operational.  If  feasibility  were  demonstrated, 
the  experiment  could  be  programmed  to  collect  additional  background  data  within 
desired  latitudinal  regions  and  at  desired  points  within  the  diurnal  and  seasonal 
cycles.  In  addition,  the  experiment  is  highly  flexible  with  respect  to  wavelength 
coverage,  so  that  data  could  be  collected  at  specified  wavelengths  to  potentially 
identify  phenomenology  that  may  be  controlling  the  fine  scale  clutter. 
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With  this  high  potential  payoff,  a  program  was  undertaken  to  assess  the 
feasibility  of  extracting  fine-scale  spatial  clutter  information  from  the  Polar 
BEAR/AIRS  Imager  data. 

The  AIRS  UV  Imager  has  been  measuring  aurora  and  airglow  data  from  a  nearly 
circular,  non-sun  synchronous,  polar,  1000  km  elevation  orbit  since  December  1986. 
The  instrument  is  mounted  in  an  earth  nadir-viewing  direction  on  the  Polar  BEAR 
Satellite,  a  Navy  Transit  class  satellite  that  is  a  3-axis  passively  stabilized 
platform  with  a  momentum  wheel. 

The  instrument  has  significant  measurement  versatility.  The  sensor  can  operate 
in  either  a  multispectral  imaging,  nadir  spectroscopy,  or  multispectral  nadir 
photometry  mode,  selectable  by  ground  command.  Regional  coverage  is  also 
selectable,  determined  by  the  positions  of  receiving  stations  and  operational 
commands.  Three  fixed  ground  receiving  stations  are  located  in  Greenland,  Norway, 
and  Fort  Churchill,  and  there  is  also  a  transportable  receiving  station  that  can  be 
positioned  where  needed  to  access  a  desired  orbit.  Measurements  are  generally  made 
north  of  45°N  due  to  power  constraints.  Low  altitude  data  can  be  collected  by 
adjusting  the  latitude  of  operation. 

It  is  the  photometer  mode  data  that  provides  the  potential  for  extracting 
fine-scale  spatial  clutter  information.  In  the  photometer  mode  the  AIRS  UV  Imager 
is  in  a  fixed  nadir  viewing  orientation.  The  sensor's  f ield-of-view  (5  km  across 
track  x  20  km  along  track  IFOV)  moves  along  the  ground  track  at  a  rate  of 
6.37  km  s“l.  Measurements  are  made  with  an  integration  period  of  6.8330  ms,  during 
which  the  IFOV  slides  -A3. 5  meters  along  the  ground  track.  Over  a  cycle  time  of 
2.9987  s  a  series  of  326  measurements  are  made,  and  these  are  summed  in  the 
conventional  processing  scheme  to  yield  resultant  integrated  5  km  x  20  km  pixel 
radiances . 

This  paper  describes  the  methodology  and  the  results  obtained  at  1304A 
wavelength  from  an  analysis  of  the  AFGL  Polar  Bear  experiment.  The  basic 
measurement  equipment  provided  data  of  a  spatial  resolution  of  20  km  over  a  large 
portion  of  the  earth.  The  instrumentation  also  provided  sampled  outputs  as  the 
footprint  scanned  along  the  measurement  track.  The  combination  of  the  fine 
scanning  and  large  area  coverage  provided  opportunity  for  a  spatial  power  spectral 
analysis  that  in  turn  provided  a  means  for  extrapolation  to  finer  spatial  scale;  a 
companion  paper  discusses  the  physical  basis  for  this  extrapolation. 

2.  RELATIONSHIP  OF  MEASURED  POLAR  BEAR  DATA  TO  SPATIAL  CLUTTER  SPECTRA 

The  Polar  Bear  Photometer  mode  data  was  collected  with  an  equivalent  ground 
footprint  of  5  x  20  km  that  was  scanned  along  the  nadir  projection  of  the  satellite 
trajectory.  The  detector  output  at  time  to  can  be  modeled  as  I(to)  where 

to 

mo)  ■  |  {  J  J  LB(x,y)  W(x-vxt,y)  dxdy]  dt 

to-AT 


where 


Lg  is  the  background  radiance  level 


AT  is  the  integration  time  for  the  measurement 
K  is  the  instrument  calibration  constant 
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and 


W(x>y)  is  the  spatial  convolution  of  the  instrument  point  spread  function  with 
the  5  x  20  km  rectangular  footprint  of  the  sensor. 

The  spatial  variability  of  the  background  radiances  Lg  is  traditionally 
characterized  by  a  spatial  Power  Spectral  Density  Sg(Wx,Wx).  One  of  the  great 
utilities  of  this  characterization  is  the  fact  that  the  mean-squared  value  of  the 
variation  of  the  measured  signal  in  any  pixel  of  size  ix  by  ly  can  be  obtained  from 
Sg  via  the  integration 


(IJJ 
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where 

e C * )  is  the  expected  value  of  (•) 


J  is  the  radiant  intensity  (watts/sr)  observed  by  a  sensor  with  pixel  size  £x 
by  ly 


Sg  is  the  power  spectral  density  in 


(w  t  c.  I  srl2 


(1  |  cru. 

wx,  Wy  are  the  spatial  frequency  variables 


and  we  have  assumed  that  the  optical  transfer  function  of  the  sensor  can  be  modeled 
by  a  Gausian  point  spread  function  with  radius  r^. 


Note  that  in  the  situation  where  fcy  -*■  o  this  relationship  can  be  approximated 
by 

t!  I  j  1 2)  =  cyy>4  J  /  sb  '"x’”,5  dux% 

or  the  variance  of  the  radiant  intensity  is  a  function  only  of  the  pixel  size. 


This  relationship  suggests  that  the  power  spectral  density  characterizing  the 
background  can  be  obtained  by  Fourier  transforming  the  measured  Polar  Bear  output. 
Since  this  data  is  in  discrete  form,  i.e.,  values  at  the  end  of  each  sample 
interval  AT,  one  utilizes  the  discrete  Fourier  transform 
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F  [ I(nAT)  ]  =  F  (^)  = 
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l 

n=l 


I(nAT) 


-J(»4T)  ) 


where 

k  is  the  discrete  Fourier  frequency  index,  and 

N  is  the  total  number  of  time  samples  transformed. 

One  can  then  obtain  the  following  relationship  between  this  transform  of  the  data 
and  Sg 
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where  the  data  calibration  factor  K  has  units  of  output  intensity  (perhaps  counts 
in  a  photomultipler)  for  a  measurement  pixel  of  dimensions  lx,ly,  per  radiance 
level  (w  |  em^  |  sr) . 

It  is  clear  that  the  Fourier  transform  of  the  measured  data  does  not  yield  the 
desired  background  power  spectral  density  directly  but  that  there  is  an  integrated 
weighting  due  to  the  nature  of  the  measurement  pixel  size  lx,  ly,  its  optical  point 
spread  function,  a  processing  function  which  depends  on  the  number  of  samples  used 
in  the  discrete  transform,  and  the  net  linear  motion  of  the  footprint  (VXAT)  that 
occurs  during  the  sampling  time  AT.  In  certain  cases  this  integrated  weighting  can 
be  deconvolved  to  obtain  the  desired  S3.  For  example,  if  Sg  can  be  written  as  a 
function  of  wx  times  a  function  of  Wy,  i.e.,  Sg(wx,Wy)  =  f(wx)  f(wy)  then  as  the 
number  of  the  Fourier  samples  is  increased  and  the  footprint  motion  decreased  so 
that  NVXAT  grows  large  then 
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Thus  the  desired  PSD  Sg  can  be  obtained  if  the  optical  parameters  of  the  sensor  are 
known,  i.e.,  the  footprint  ix  and  blur  circle  radius  r^. 

A  potential  equally  difficult  problem  involved  with  the  extraction  of  the 
background  PSD  from  the  measured  data  is  the  noise  associated  in  the  sensor.  In 
the  Polar  Bear  analysis  being  discussed  the  limiting  noise  is  the  shot  noise  or 
counting  noise  associated  with  the  limited  number  of  optical  photons  available  in 
the  receiver  during  the  sample  time  AT.  The  contribution  of  the  shot  noise  to  the 
Fourier  transform  of  the  measured  data  can  be  assessed  by  assuming  that  the 
probability  of  the  number  of  counts  n^  in  the  i^*1  sample  period  can  be  described  by 
a  Poisson  distribution  with  probability  density 

-h  Mk 

Prob  (n.  =  k}  =  e  “ 

where  M  is  the  average  number  of  counts  in  the  sample  period,  note  that  the  RMS 
value  of  the  counts  in  then  /M.  The  expected  value  of  the  Fourier  transform  of  a 
sequence  of  N  such  samples  satisfying  this  statistical  law  is  then  determined  by 

e|  1  F[n.]  |  2}  =  N2M2  +  KM  at  zero  frequency 
NM  at  all  other  frequencies 

Analysis  of  a  sequence  of  approximately  3  seconds  of  Polar  Bear  data  (with  M  =  256) 
resulted  in  a  normalized  Fourier  transform  spectrum  shown  in  Figure  1.  The  plot  is 
of  the 


SPATIAL  FREQUENCY  (CY  /  M) 

Figure  1.  Spectral  Analysis  of  a  Three  Second  Segment  of  Polar  Bear 
Photometer  Data 
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ZlHjamii  2 

N2(Ave  I)2 

i.e.,  the  square  of  the  transform  divided  by  the  square  of  the  product  of  N  and  the 
average  value  of  the  measured  data.  On  the  same  curve,  the  normalized  level  of  the 
Fourier  spectrum  of  a  purely  shot  noise  process  is  also  shown.  It  is  clear  that 
the  measured  data  is  dominated  by  the  shot  noise  in  the  observation.  The 
likelihood  of  extracting  any  useful  information  from  such  a  sample  is  very  small. 
Fortunately  as  larger  sequences  of  data  were  analyzed  one  does  obtain  useful 
information  at  lower  spatial  frequency,  for  example  Figure  2  shows  the  Fourier 
spectra  of  a  sequence  of  N  —  32768  points  (approximately  300  seconds  of  data)  which 
shows  the  shot  noise  level  and  the  component  of  the  measurement  due  to  the  spatial 
variation  of  the  background  that  emerges  at  lower  frequencies. 

3 .  Extrapolation  of  Low  Frequency  Spectral  Data 


The  procedure  that  has  been  evolved  to  extract  the  background  radiance  PSD,  Sg, 
from  the  Fourier  spectrum  of  the  measured  data  is  based  on  an  assumed  form  for  Sg 
namely  that 
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One  can  then  demonstrate  that  if  the  discrete  Fourier  transform  variable  k  is  such 
that 
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then  the  fourier  spectrum  of  the  measured  data  can  be  approximated  as 

2  NV  AT  •• 
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Figure  2.  Spectral  Analysis  of  Three  Hundred  Second  Segment  of  Polar  Bear 
Photometer  Data 

A  A 

Moreover,  in  a  footprint  lx  such  that  lx  «  L,  the  variance  of  the  measured  data 
can  be  approximated  by 


E[lunlT)-?|2!=^ 

2tt  L 

A 

wher^  K  is  the  instrument  calibration 
K  =  lx/lx  K 


6+1 

2  (6-1) 

A 

constant  for  a  footprint  dx;  note  that 


These  two  expressions  suggest  the  following  procedure  for  obtaining  the  desired 
constants  (6,  A,  and  L) ,  defining  the  background  variance  Sg,  from  the  measured 
spectra 


1)  in  a  region  of  the  measurement  spectra  not  dominated  by  shot  noise  establish 
a  best  fit  of  the  function 


-  _  constant 

Thus  establishing  6  and  one  relationship  between  A  and  L  given  by 
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where  |  F0  j  2  is  some  point  on  the  spectrum  corresponding  to  the  variable  kQ  where 
the  spectra  and  the  fitting  function  coincide  (assumed  to  be  at  as  large  a  value 
of  kQ  as  possible  but  below  the  region  dominated  by  shot  noise). 


2)  utilize  the  approximate  expression  for  the  variance  to  establish  a  second 
relationship  between  A  and  L.  These  two  relationships  can  be  combined  to 
obtain  expressions  for  L  and  A  given  by 


VT  k  &  [L  (1  +  z2)B+1/2  W  [  I  F  I  2  -  NI] 


and 


x  o 


'A  -  [cf  |f(„4T,  -  T-l  2J  -?1 

K  (B+l) 


A  ^ 

where  I  is  the  equivalent  measured  data  in  a  footprint  £x  corresponding  to 


N  V  AT 

k  =  - *— 

°  l 

x 

Note  that  the  measured  Polar  Bear  data  was  repeatedly  sampled  at  time  steps 
corresponding  to  increment? 1  motion  of  the  sensor  footprint  (VxAt).  Thus  if  one 
wanted  an  £x  corresponding  to  "M”  multiple^  of  the  basic  Polar  Bear  footprint  (£x) 
then  samples  of  data  corresponding  to  the  £x  footprint  can  be  obtained  by  the 
indexed  sum 

I  =  I (n)  +  I(n  +  k)  =  I(n  +  2k)  +  -  +  I[n  +  (M  -  l)k] 


of  the  basic  Polar  output  data  which  has  k  data  points  within  each  footprint,  i.e. 
k  VXAT  =  £x.  Note  also  that  these  expressions  have  been  corrected  for  the  shot 
noise  contribution  to  both  the  Fourier  spectra  of  the  measured  data  and  the 
variance  of  the  measured  data.  In  addition  the^constant  K  must  be  the  equivalent 
instrument  scale  factor  for  the  footprint  size  ix  (ft  =  MK  for  the  above  situation). 


These  results  can  also  be  applied  directly  to  obtain  a  measure  of  the  RMS 
level  of  the  background  clutter  since 


°i  l  =  ra<liant  intensity  =  £^£ 

x  y  X  ^ 


if  £x  and  ly  are  both  much  less  than  the  model  constant  L.  In  particular  the  RMS 
level  of  the  background  radiance  is  given  by 
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4 .  Example  Results  of  the  Polar  Bear  Data  Analysis 


Using  the  procedure  developed  in  Section  3,  a  series  of  Polar  Bear  photometer 
mode  data  measured  at  1304  A  was  analyzed  to  obtain  estimates  of  the  background 
clutter  power  spectral  density.  A  typical  curve  fitting  to  the  measured  spectra 
obtained  for  N  =  32,768  data  pounds  or  approximately  300  seconds  is  shown  in 
Figure  3  from  which  one  obtains  the  slope 

8  s  1.97 


and  the  fitted  point  kc  =  24  (corresponding  to  a  spatial  frequency  of  1.65  x  lO-3 
cycles/mcter  and  a  footprint  size  £x  of  60  km)  for  |  F0  |  2  =  1.316  x  107.  For  the 
same  sequence  of  data  =  0.329  x  1033.  The  instrument  calibration  constant  K 
was  determined  to  be 


6.06  x  1010 


counts 
.  2 

w/cm  -sr 


The  variance  of  the  measured  data  in  the  60  km  equivalent  footprint,  i.e., 
j  1  -  ?  |  2}  was  obtained  as  (14.5)2  counts2  and  the  mean  number  of  counts 
I  =  1.66.  The  value  of  L  obtained  from  the  corresponding  expressions  in  Section  3 
was  2  x  103  km  which  for  the  sake  of  physical  reference  is  approximately  one 
twentieth  the  circumference  of  the  earth. 
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Figure  3.  Example  of  the  Data  Fitting  Procedure  for  Polar  Bear  Data 
(Summer,  Mid  Day  Observations  at  300  Seconds  at  Latitude  60  -  75°  North) 
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Similar  analysis  was  conducted  at  different  latitudes  with  results  shown  in 
Figure  4  for  the  slope  0,  the  RMS  value  of  the  spectral  radiance,  and  the  ratio  of 
the  RMS  to  mean  radiance.  Note  that  the  nominal  spectral  half  width  of  the 
measured  instrument  was  35  A  which  allowed  the  results  to  be  quoted  as  spectral 
radiance. 

The  resulting  background  radiance  levels  and  the  corresponding  relationship  for 
the  RMS  level  of  the  clutter  in  the  footprint  (a  -(constant)  £x£y)  was  justified 
for  the  extrapolation  procedure  applied  to  this  data,  i.e.  the  approximations 
employed  in  the  analysis  are  met  by  the  data.  Thus  one  has  a  means  for  estimating 
clutter  at  finer  spatial  scale  than  the  Polar  Bear  footprint.  However,  since  the 
measured  data  at  finer  scales  is  contaminated  by  counting  noise  one  can  not  be  sure 
that  the  extrapolation  procedure  to  scales  such  as  10-100  meters  is  valid,  i.e., 
there  may  be  physical  processes  that  enter  at  the  scale  which  invalidate  the 
fractal  nature  of  the  phenomena.  The  discussion  of  the  companion  paper  provides 
some  justification  for  the  extrapolation  into  their  region  based  upon 
phenomenological  considerations.  It  is  clear,  however,  that  direct  measurements  at 
these  scale  lengths  are  necessary  in  order  to  establish  the  validity  of  these 
cone lus ions . 
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Figure  4.  Polar  Bear  Measured  Spatial  Radiance  Clutter  Levels  of  Sunlit  Earth 
Background  at  1304  A 


